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Abstract 

Fourier transform of the generalized parton distributions (GPDs) at zero skewness with respect 
to the transverse momentum transfer gives the distribution of partons in the impact parameter 
space. We investigate the GPDs as well as the impact parameter dependent parton distributions 
(ipdpdfs) by expressing them in terms of overlaps of light front wave functions (LFWFs) and 
present a comparative study using three different model LFWFs. 
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I. INTRODUCTION 



Deeply virtual Compton scattering (DVCS) 7*(g) + p{P) l{,<i') +p{P'), where the 
virtuality of the initial photon = —q^ is much large compared to the squared momentum 
transfer t = —{p — p'Y , provides a valuable probe to the structure of the proton near the 
light cone. At leading twist, QCD factorization holds and DVCS amplitude can be expressed 
as a convolution in x of the hard 'y*q —>■ "fq Compton amplitude with the generalized parton 
distributions (GPDs) Here x is the light cone momentum fraction of the active quark. 
The skewness ( = measures the longitudinal momentum transfer in the process. 

GPDs are richer in content about the hadron structure than ordinary parton distributions 
(pdfs). On one hand, x moment of the GPDs give hadron form factors measurable in 
exclusive ep —>■ ep scattering, on the other hand in the forward limit, i. e. for zero momentum 
transfer, they reduce to ordinary pdfs measurable in inclusive processes. Thus they provide 
a unified picture of the hadron. 

GPDs are off-diagonal overlaps of light-front bilocal operators and unlike ordinary pdfs, 
they do not have an interpretation of probability densities. It has been shown in [2] that 
a Fourier transform (FT) of the GPDs with respect to the transverse momentum transfer 
A_|_ at zero skewness ( gives the distribution of partons in the transverse position or im- 
pact parameter space. They are called impact parameter dependent parton distributions 
(ipdpdfs) q(x, b±). The impact representation of pdfs on the light front was first introduced 
by Soper jsj in the context of the FT of the elastic form factor. q{x, b±) give simultaneous 
information on the distribution of quarks as a function of x and the transverse distance b± 
of the parton from the center of the proton in the transverse plane. Ipdpdfs obey certain 
positivity constraints and thus, it is legitimate to physically interpret them as probability 
densities. In fact, this interpretation is not limited by relativistic effects in the infinite mo- 
mentum frame. q{x, b±) are defined for a hadron state at sharp momentum P"*" localized 
in the transverse plane such that the transverse center of momentum is at R± = (one 
can also work with a wave packet state localized in the transverse position space in order 
to avoid the state to be normalized to a delta function). When the target is transversely 
polarized, the distribution of partons in the impact space is no longer axially symmetric 
and the deformation is described by the FT of the GPD E{x, 0, t). This distortion has been 
shown to be connected with Sivers effect 0, S]. A 3-D picture of the quarks and gluons 
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in the proton has been proposed in in terms of reduced Wigner distributions which are 
related to the FT of the GPDs in the rest frame of the proton. In [7] it has been shown 
that a FT in the skewness ( provides a boost invariant longitudinal position space picture of 
the proton and is analogous to optical diffraction obtained in single slit experiment. Thus, 
GPDs provide a unique picture of the hadron in transverse and longitudinal position space. 

GPDs can be expressed in the light-front gauge as overlaps of the light-front wave func- 
tions (LFWFs) of the target hadron. These are off-forward overlaps in general and one 
requires not only a particle number conserving n — > n overlap similar to forward pdfs but 
also ann-|-l— s>n — 1 overlap, when the parton number is decreased by two. When ( is 
zero, the second contribution vanishes. GPDs and the ipdpdfs have been investigated in 
several phenomenological models, for example in chiral quark model for the pion in the 
constituent quark model [ol, in terms of a power law ansatz of the light cone wave function 
for the pion (lo|, in the context of investigating the color transparency phenomena in 
the transverse lattice framework for the pion 121] and in the lattice framework 13|. 

In this work we present a comparative study of the GPDs as well as the ipdpdfs using 
using several phenomenological models of hadron LFWFs. In [l^ DVCS amplitude at one 
loop in QED has been computed for a fermion. In effect, one represents a spin 1/2 system 
as a composite of a spin 1/2 fermion and a spin-1 vector boson, with arbitrary masses 15 1. 



Similar models have been used in 



16 



171]. This one loop model is self consistent since it has 



the correct correlation of different Fock components of the state as given by the light-front 
eigenvalue equation. In [7] a simulated model for the meson has been derived from the above 
LFWF, by taking a derivative with respect to {wrt) the bound state mass square M^, which 
appears in the denominator of the wave function, thus improving the behaviour of it near the 
end points at a; = 0, 1. In this model, DVCS amplitude is purely real. A similar power law 
behaviour has been used in [l^ to construct the GPDs for a meson. Differentiating once 
wrt generates a mesonlike k± behaviour and differentiating wrt internal fermion mass 
and the arbitrary gauge boson mass simulates a protonlike behaviour of the LFWFs. 
We present results in both these models. Models for LFWFs of hadrons in 3 -|- 1 dimensions 
displaying confinement at large distances and conformal symmetry at short distances have 
been obtained using ADS/CFT method. We also present ipdpdfs in this model. 
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II. GENERALIZED PARTON DISTRIBUTIONS 



The kinematics of DVCS process is given in detail in Q . We work in the frame of 
We take skewness ( to be zero. Momentum transfer is purely transverse; 
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(2.1) 



The generalized parton distributions H, E are defined through matrix elements of the bi- 
linear vector currents on the light-cone: 

dy- 



y+=0,yj_=0 



1 



2P^ 



UiP') 



2M 



UiP) 



(2.2) 



here P = |(P' + P) is the average momentum of the initial and final hadron. 

The off-forward matrix elements can be expressed as overlaps of the light front wave 
functions For non-zero skewness ( there are diagonal parton number conserving con- 
tributions in the kinematical region ( < x < 1 and ( — 1 < x < 0. There are off diagonal 
parton number changing contributions in the region < x < (. In our case C = and 
the only relevant kinematical region is < x < 1. These correspond to the target helicity 
non-flip (++) and helicity flip (H — ) contributions, respectively. If we consider a spin 1/2 
target state consisting of a spin 1 particle and a spin 1/2 particle, these contributions can 



be expressed in terms of the 2-particle LFWFs 
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where 
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2M 
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l-x) A, . 



(2.4) 



(2.5) 



Here '?/'aiA2(-^' ^1) lowest (two-particle) Fock component of the hadron LFWF with 

helicity up and Aj, z = 1, 2 are the intrinsic helicities of the internal particles. 
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FIG. 1: Generalized parton distributions in model 1; (a) H{x,0,t) and (b) E{x,0,t) as a function 
of X for fixed values of in units of MeV. 



The impact parameter dependent parton distributions are defined from the GPDs by 
taking a FT in A^, 

1 



q{x,b±) 



-ib I -A 
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(2.6) 



where b± is the impact parameter conjugate to A_ 



III. SIMULATED MODEL CALCULATIONS 



In this section, we calculate the GPDs in simulated mode 



with the two-particle wave function for spin-up electron 
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s o: 
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hadron LFWFs. We start 
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FIG. 2: Generalized parton distributions in model 2; (a) H{x,0,t) and (b) E(x,0,t) as a function 
of X for fixed values of in units of MeV. 



Following the same references, we work in a generalized form of QED by assigning a mass 
M to the external electrons and a different mass m to the internal electron lines and a mass 
A to the internal photon lines. The idea behind this is to model the structure of a composite 
fermion state with mass M by a fermion and a vector 'diquark' constituent with respective 
masses m and A. Similarly, the wave function for an electron with negative helicity can also 
be obtained. In Eq. (13.21) . the bound state mass M appears in the energy denominator. 
As discussed in j?], a differentiation of the QED LFWFs with respect to improves the 
convergence of the wave functions at the end points: a; = 0, 1, as well as improves the kj_ 
behaviour, thus simulating a bound state valence wavefunction. Differentiating once with 
respect to will generate a meson-like behaviour of the LFWF. However, this is not a 
model for a meson wavefunction since the two constituents have spin half plus spin one. If 
we differentiate once more we simulate the fall-off at short distances which matches the fall- 
off wavefunction of a baryon, in the sense that the form factor Fi{Q^) computed from the 
Drell-Yan- West formula will fall-off like ^ . Here we have the analog of a two-parton quark 
plus spin-one diquark model of a baryon, not three quarks. Overlaps of these wavefunctions 
in the same way as have been done for the dressed electron wavefunctions will simulate the 
corresponding GPDs. This is the approach we are following here. 

We differentiate (p{x, k±) with respect to the bound state mass in order to simulate 



6 




1 1 1 1 1 

(b) 


b=0.05 




b=0.1 






--- b=0.5 




— ^ 'I , 1 , 1 


^ T — _ 



0.2 



0.4 



0.6 



0.8 



FIG. 3: Impact parameter dependent parton distributions in model 1; (a) q{x, b±) and (b) e(x, b±) 
as a function of x for fixed values of I 6_l I in units of MeV~^ 



the LFWFs of a meson-like hadron. In other words, we take 

\ X 1—x J 

where N is the normalization constant. We normalize the 2- particle LFWF to 1. 

The GPD in this model when the target helicity is not flipped (helicity non-flip) becomes 
(from Eq. flOD ): 

r/ r 

H{x,t) = ^ArV(l-a;)[((l + a;2)/i + (l+x')/2 + {(l + a;')P 

+ 2{l-x)\Mx-myh])], (3.4) 



where 
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LiLi Jo ^ 
3 = I TTF2=^^l dp ^3 , (3.5) 



here Li = m2(l - x) + A^x - M'^x{l - x), L2 = k\- 2k^ ■ A_l(1 - x) - B and 

D = f3m\l -x)- f3M^x{l - x) + f3X^x - (1 - f3)B - (1 - (3)\l - xfAl, (3.6) 



7 




FIG. 4: Impact parameter dependent parton distributions in model 2; (a) q{x, b±) and (b) e(x, b±) 
as a function of x for fixed values of | 6_l | in units of MeV^^. 

where B = M^a;(l— x)— A5_(l— x)^— m^(l— x)— xA^ and P = i?+M^a;(l— x)— m^(l— x)— A^x. 
The hehcity flip part as obtained from Eq. 12.41 is given as: 

E(x, t) = 2MN^—x^{l - xf{Mx - m)h. (3.7) 

Next, we simulate the LFWF for a proton-like hadron by differentiating eq. (13. 2p first 

wrt m? and then wrt A^. We get, 

r , , (9¥^(x,fc^) 2e AT 

^ ^^'^^ ' 5mW ' x(l-x)3/2 /^2_a±!!!!_^i±^^' ■ ^ ^ ^ 

Here N" is the normalization constant. As before, we normalize the two particle LFWF to 
1. 

The helicity non-flip GPD in this model becomes : 

H{x,t) = — ArV(l-x)[((l + x')j; + (l + x2)/^ + {(l + x')P 

+ 2(1 -x)2(Mx -m)2/3})], (3.9) 



where 



= Jm^ Jo ^^^D^ 



LlLl Jo ' 
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FIG. 5: Impact parameter dependent parton distributions q{x, b±) in (a) model 1 and (b) in model 
2 as a function of | 6^ | in units of MeV~^ for fixed values of x. 

The helicity flip part is given by, 

E{x, t) = 2MAf^—x%l - xf{Mx - m)I'^. (3.11) 

Here we refer to the first model LFWF as model 1 and the second as model 2. In both 
cases, we take M = 150 MeV, m = X = 300 MeV [t]. In fig. 1 (a) we have plotted 
the helicity non flip GPD H{x, 0, t) as a function of x for fixed values of t. An interesting 
aspect of these models is that in both of them one can generate a helicity flip contribution 
E{x,0,t), although the k± behaviour of model 1 is similar to a meson. This is because, 
we have simulated this k± behaviour by using a two-particle composite state, where one 
of the constituents has spin 1/2 and the other has spin 1. We have plotted the helicity 
flip GPD E{x,0,t) in model 1 in flg. 1 (b). In flgs 2 (a) and (b) respectively, we have 
shown H{x,0,t) and E{x,0,t) for model 2. From the flgs, one can see that the qualitative 
behaviour of the GPDs in both models are similar. The GPD H{x,0,t) increases with x, 
reaches a maximum and then falls to zero at a; 1 independent of the momentum transfer 
t. As X is the momentum fraction of the active quark and at x = 1 the active quark carries 
all the momentum, contributions from the other partons are expected to be zero in this 
limit and H{x,0,t) is expected to become t independent. The peak of H occurs at higher 
values of x, and it shifts towards even higher values of x as | t | increases, which means 
that the active quark is more likely to have a larger momentum fraction. In fact, here we 
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FIG. 6: Impact parameter dependent parton distribution in the holographic model. LHS is for the 
ground state L = 0, k = I and RHS is for the first exited state L = 1, k = 1. Aqcd = 0.32 GeV 
and b± is given in GeV~^. 

are considering only the leading Fock space component of the target hadron state. Higher 
Fock space components are more likely to contribute in the small x region. However, the 
magnitude of the peak decreases for high \ t \, which is expected as the form factor decreases 
with I t I for high t. In model 2, the peak shifts to larger x faster with an increase of | t |. The 
GPD E{x,0,t) is associated with the amplitude when the nucleon helicity flips and quark 
helicity does not flip. This is known to be related to the quark orbital angular momentum. 
The functional dependence of E{x, 0, t) is different from H{x, 0, t) in both models. The peak 
is at X = 0.5 for smaller values of | t | but shifts to larger values of x as | t | increases. It is 
zero at both a; = 0, 1, same as H{x, 0, t). Like H{x, 0, t), it also becomes independent of | t | 
as X ^ 1. E{x,0,t) decreases as | t | increases, basically as the first moment of E{x,0,t) 
gives the form factor F2{t), which decreases for large \t\. 

Taking a Fourier transform (FT) of the GPDs with respect to the transverse momentum 
transfer A^, one gets the ipdpdfs. In figs 3 and 4, we have plotted the ipdpdfs q{x, b±) and 
e{x, b±) defined as in eq. (12. 6p for models 1 and 2 respectively. The ipdpdfs show the same 
qualitative behaviour in both models, which is expected as the GPDs themselves behave in a 
similar manner. In other words, qualitatively, parton distributions in the impact parameter 
space are not that sensitive to the large k± behaviour of the LFWFs. The peak is shifted 
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to larger x in model 1. q{x,b±) goes to zero at x = and x = 1. In fact, in the limit 
X —>■ 1 the active quark with momentum fraction x is very close to the transverse center of 
momentum, and the transverse width of the ipdpdf vanishes [2]. The peak in x decreases as 
bj_ increases, thus the transverse profile in bj_ is peaked at | 6_l |= and falls away further 
from it. The ipdpdf with the nucleon helicity flip, e(x, b±) is especially interesting, as for 
a transversely polarized target, it measures the distortion of the parton distribution in the 
transverse plane. In addition, the second moment of e{x, b±) gives how the distribution of 
momentum of quarks of a particular flavour in the transverse plane changes when the nucleon 
is not polarized longitudinally. This receives contribution from quarks with nonzero orbital 
angular momentum and is shown to be connected with Sivers effect. e{x,b±) is peaked at 
X = 0.5 in both models and the height of the peak decreases as | 6^ | increases. In Fig. 5, 
we have shown in ipdpdfs q{x, b±) in both models as a function of | 6_l | for fixed values of 
X. It is interesting to see that in both models, there is a primary peak at | 6^ | = and then 
there are several secondary peaks. Positions of these secondary maxima are independent of 
X. It is to be noted that in model 2, the curve for x = 0.5 lies above the one for x = 0.8, 
whereas in model 1, the situation is opposite. This is because in model 1, the peak is shifted 
towards higher values of x than in model 2 (see figs 3 and 4). 



In 



20], it has been shown that the string amplitude $(-2) defined on the fifth dimension 



in ADS^ space can be mapped to the light-front wave functions of hadrons in physical 
spacetime. The holographic variable z corresponds to the impact variable ^ where = x(l — 
x)b']_, b± being the impact parameter. The effective 4-dimensional Schroedinger equation for 
the bound states of massless quarks and gluons exactly reproduce the ADS/CFT results and 
gives a realistic description of the light quark meson and baryon spectrum. The normalized 
two particle LFWF for the bound state is given by. 



*L,fc(x,6±) = BL,k^xil-x)JLi^(3L,kAQCD); (3.12) 

where BL,k = ^QCD[i-^)^TrJi+LiPL,k)Ji-Lil3L,k)]~^^^, (3L,k is the fc-th zero of Bessel function 
Jl and Aqcd = 0.32 GeV. For ground state L = 0, k = 1 and we have 

^oAx,b^) = Kqcd4^^)^-^^^B^^. (3.13) 

In fig. 6 we have plotted the ipdpdf for the ground state as well as the first exited state 
as functions of x and | 6± |. One has to note that in this model b^ cannot be very large 
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as 2; < j^^- The primary peak of the ipdpdfs occur for the ground state (L = 0, = 1) 
at I 6_L 1= 0. In the first exited state (L = 1, A; = 1) it is shifted away from | &_l |= 0. In 
both cases, there is a secondary peak at a higher value of | 6_l | for certain x values. Similar 
secondary peaks have been observed for the other two models considered, as we already 
noted. 

IV. SUMMARY AND CONCLUSION 

We have investigated GPDs as well as impact parameter dependent parton distributions 
at C = in three different phenomenological models of hadron LFWF. We simulated models 
for meson-like and proton-like hadrons respectively, starting from the two-particle LFWF 
of a dressed electron in a generalized form of QED, and differentiating the denominator 
wrt the squared masses. This improves the behaviour at the end points x = 0, 1 as well as 
changes the k± behaviour. We found that the GPDs and the ipdpdfs are not that sensitive 
to the k_i_ behaviour of the LFWFs. The peak of the ipdpdfs in x for fixed b± is shifted 
towards higher x in the simulated model for a meson-like hadron. We have also investigated 
the ipdpdfs in the recently proposed QCD holographic model, and shown the behaviour over 
a larger 6^ range. In all three models, ipdpdfs show several secondary peaks in b± besides 
the primary one. 
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